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Abstract. Using photometric selection criteria, and in a final step visual inspection of spectra, we find 65 DQ 
white dwarfs in the Data Release 4 of the Sloan Digital Sky Survey. 37 of these are new identifications, three 
more are reclassified DQ from DC, the other 25 have been described as DQ before in the literature. We derive 
effective temperatu res and carbon abundances for 60 of these stars. The results confirm the main conclusions of 
iDufour et alJ <l2005f) . The majority of stars defines a clear sequence in the C abundance vs. T e g diagram, with high 
abundances found at high T e g and vice versa. We also confirm with high significance a second sequence with an 
abundance about f dex higher at the same T e g, and discuss the nature and possible origin of the high-C sequence. 
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1. Introduction 

White dwarfs of spectral type DQ are defined as show- 
ing atomic or molecular features of carbon in a he- 
lium atmosphere. Contrary to other metals (notably Ca) 
found in cool white dwarfs, which are thought to be 
provided through an accretion process from outside the 
star, the most widely accepted explanation for the DQs 
is dredge-up from the underlying carbo n/oxygen core 
through the expanding He convection zone l|Koester et all 
ll982t|Pelletier et allll98ri . A detailed summary of the ba- 
sic properties of this class and earlier references to spectro- 
scopic analyse s has been given rec ently in a comprehen- 
sive paper by IDufour et al ] l)2005l = DB05 henceforth) 
and need not be repeated here. 

In addition to the general understanding of the nature, 
evolutionary relations and origin of the peculiar spectral 
types of white dwarfs, the DQs are of special interest, 
because they provide information about the deeper lay- 
ers of the stars. The transition between the outer he- 
lium layer and the carbon/oxygen core - the result of 
the helium burning in the progenitor - is not abrupt, but 
rather gradual as det ermined by the diffusion equilibrium 
l|Koester et allll982j) . We can predict the struc ture of this 
transition quantitatively ijPelletier et all Il986j) , as well as 
the depth of the helium convection zone at a given ef- 
fective temperature and mass of the star. The observed 
carbon abundance in the atmosphere is then directly re- 
lated to the thickness of the helium layer remaining from 
the previous nuclear evolution, which is a very important 
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parameter for the understanding of stellar evolution. The 
first estimate of this thickness by IPelletier et all l|l986^ 
- logMiie/M w —3.75 - was much thinner than pre- 
dicted by evolutionary calculations of the AGB phase. 
However, the recent abundance determinations by DB05, 
together with the n e w mo dels of the outer layers by 
iFontaine fc Brassard! l)2005|) are in much better agree- 
ment with theoretical expectations. Nevertheless, there 
are still many problems understanding the detailed con- 
nection of the DQ stars to possible he l ium-r ich progenitors 
l|Althaus et all l2005t IScoccola et al"ll2006|h 

When searching for correlations of DQ parameters, 
e.g. between atmospheric carbon abundance and effec- 
tive temperature, it is very helpful to have available a 
homogeneous sample of observations. It has been found 
in previous studies that abundances determined from 
atomic lines in the UV or molecular ba nds in the optical 
may sometimes differ significantly fe.g. IProvencal et all 
2002j), possibly blurring any correlation if only one or 
the other observation i s ava ilable for different objects 
( Weidemann fc Koesteii Il995lh Such a homogenous sam- 
ple is currently provided by the Sloan Digital Sky Survey, 
and DB05 provides the first quantitative analysis of 40 
DQs from the Data Release 1 (DRI). In this paper we 
report similar results for 40 n ew DQs (and 20 previously 
know n) extracted from DR4 l|Adelman-McCarthv et all 
l200fi|h 



2. Selection of DQ candidates 

DR4 of the SDSS contains spectra for approximately 
240000 objects classified as "point source" or "unknown" . 
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In order to narrow down this huge database to tractable 
numbers we have applied three consecutive selection steps, 
the first two based on the SDSS photometry, the last on 
the spectra. 

2.1. DQ regions in SDSS color space and theoretical 
models 

Theoretical SDSS colors can be calculated from synthetic 
spectra by convolving them with the 5 ugriz SDSS band 
passes. Our model grid is similar to that described in 
ICarollo et al"lll20f)3h . but extended to cover effective tem- 
peratures from 13000 K to 4400 K. Logarithmic C abun- 
dance by number relative to He ( = [C/He])range from —8 
to —4 in steps of 0.5. The standard grid used here has a 
surface gravity of logg = 8.00, but additional grids were 
calculated with 7.5 and 8.5 for comparison. The models 
are fully blanketed, fully including the line and molecu- 
lar band absorption in the calculation of the atmospheric 
structure. Non-ideal effects in the form of a lowering of 
the dissociation energy of the C2 molecule is included in 
a very approximate way. 

The transmission curves for the five filters (for airmass 
1.3 as recommended) were obtained from the SDSS web 
sites; the zero-points for the magnitudes were set to cor- 
respo nd to the AB system of magnitudes l|Fukugita et all 
Il996h except for small changes to u, i, z (-0.04, 0.01, 0.02) 
similar to those discussed in Eisenstein et al. (2006). 

The dependence on T e g and [C/He] is strongest in the 
u — g vs. g — r two-color diagram and much less in the re- 
maining independent indices used, r — i and i — z. Figure^ 
shows this diagr am with the theoret ical grid, and 40 ob- 
served DQs from lHarris et all 1 2003h and lKleinman et alJ 
(2004), which indeed fall in or very near the theoretical 
grid. Note that we have not applied any correction for 
galactic extinction (see DB05). 

Using this information about the position of DQs in 
the SDSS color space, we have as a first selection step 
extracted from the SDSS database all objects in the "point 
source" and "unknown" categories falling into the region 
defined by 



-0.8 <u-g< 1.5 



and 



-0.8 < g-r < 1.5, 



resulting in 86856 selected objects. This sample still con- 
tains large numbers of DA white dwarfs and quasars. We 
therefore have in a second step refined the color region us- 
ing a polygon in the u — g, g — r plane. The polygon area 
follows rather closely the outline of the theoretical grid, 
except at the hot end. Since the grid converges towards 
a single line (the models are much less dependent on the 
C abundance) the polygon at the hot end is significantly 
wider in u — g. Applying this selection step narrowed the 
candidate sample down to 44928 objects. 

2.2. Selection in 4-dimensional color space 

Since we consider for the moment only the dependence 
of colors on T e g and [C/He], it is obvious that the the- 
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Fig. 1. Theoretical SDSS color grid for DQ white dwarfs. 
Thick continuous lines are lines of constant carbon abun- 
dance from [C/He] = -8.0 (top) to —4 (bottom), step 1.0. 
Thin dashed: lines of constant T c g from 4400 K (right) to 
13000 K (left) i n steps of 200 K. Cro sses s how the observed 
DQ st ars from lHarris et all ((2003) and iKleinman et alJ 
(2004) 



oretical grid must define a 2-dimensional surface in the 
4-dimensional color space (u — g,g — r,r — i,i — z). That 
the variation is strongest in the first two components sug- 
gest that this surface might be approximated by a plane 
and this is confirmed by looking at different projections of 
the grid in two- or three-dimensional subspaces. We have 
therefore fitted a plane to the nodes of this grid using 
a x 2 minimization technique. The fit is extremely good, 
with typical deviations of grid points from the plane of 
0.03 mag. 

The equation of the plane is 
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where we have used the first two colors directly as the two 
variables for the plane. The reasonable region for DQs is 
given by the limits 



-0.55 <u-g < 1.3 



and 



0.3 < g - r < 1.4. 



This plane representation of theoretical DQ colors was 
used in the final photometric selection step. We defined a 
minimum "distance" d from the plane for each observed 
object by measuring each coordinate in units of their in- 
dividual measurement errors 
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Table 1. 40 new DQ white dwarfs. The first columns give the SDSS name and internal designation. Next is effective 
temperature determined from a fit to the photometry. The carbon abundance [C/He] is determined from a model fit 
to the spectrum, keeping the temperature fixed. The final colu mn gives a classif i cation and remarks. WD 1105+412 
and WD 0913+103 are reclassified from a DC classification in iMcCook fc SiorJ l)l999j) . GD311 is classified DC in 
iKawka fc Vennesl l)2006(l . but we believe very weak Swan bands to be present. 
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OpCCtl 3,1 1 ypC j It, CI 


SDSS .1074204.79+434835.7 


53052 


1736 


139 


7738 


-5.47 


0.04 


DQ 


SDSS .1084131.55+332915.6 


52642 


933 


16 


6810 


-6.62 


0.05 


DQ 


SDSS .1085239.66+042804.5 


52670 


1190 


177 


9555 






DQ 


SDSS J085506. 62+063904. 7 


52668 


1189 


535 


7337 


-5.99 


0.01 


DQ 


SDSS J085709. 01+060357. 4 


52668 


1189 


27 


8210 


-5.00 


0.02 


DQ 


SDSS J090449. 73+395416. 4 


52703 


1199 


595 


7319 


-5.74 


0.05 


DQ 


SDSS J090514. 78+090426. 5 


52973 


1300 


420 


8861 


-4.89 


0.01 


DQ (Ha ?) 


SDSS J090632. 17+470235.8 


52606 


898 


565 


5070 


-4.31 


0.01 


DQ (weak bands) 


SDSS J091602. 73+101110. 5 


53050 


1739 


592 


8715 


-4.81 


0.01 


DQ 


SDSS J091830. 27+484323.0 


52637 


900 


429 


8884 


-3.72 


0.06 


DQ 


SDSS J092153. 46+342136. 9 


52995 


1274 


216 


8202 


-5.46 


0.09 


DQ 


SDSS ,1092613.46+472521.1 


52637 


900 


41 


7261 


-6.41 


0.06 


DQ 


SDSS J092909. 03+331011. 7 


52991 


1593 


94 


6361 


-5.57 


0.01 


DQ 


SDSS J094014. 65+090641. 8 


52993 


1304 


45 


6169 


-7.34 


0.05 


DQpcc (rounded bands) 


SDSS J094115. 18+090154.4 


52993 


1304 


9 


9122 


-4.73 


0.02 


DQ 


SDSS ,1094138.08+441458.2 


52672 


1202 


58 


8113 


-5.42 


0.06 


DQ 


SDSS J095934. 95+453725. 4 


52703 


942 


446 


7211 


-5.52 


0.06 


DQ 


SDSS J100059. 82+100531. 7 


53053 


1308 


307 


7958 


-4.83 


0.02 


DQ, WD 1105+412 


SDSS J101750. 38+373637. 5 


52996 


1427 


216 


7497 


-5.79 


0.03 


DQ 


SDSS J101800. 00+083820. 3 


52762 


1237 


621 


7784 


-5.85 


0.01 


DQ 


SDSS J102635. 81+580714. 8 


52316 


559 


6 


8879 


-4.56 


0.03 


DQ 


SDSS J110759. 46+405910. 9 


53046 


1437 


512 


7169 


-6.45 


0.02 


DQ 


SDSS J110912. 21+424956.0 


53053 


1363 


37 


9402 


-4.84 


0.08 


DQ, WD 0913+103 


SDSS J112604. 29+441938. 6 


53062 


1365 


564 


7097 


-6.38 


0.05 


DQ 


SDSS J113534. 61+572451. 7 


53033 


1310 


485 


7385 


-6.42 


0.02 


GD 311. Kawka & Vennes f2006) 


SDSS J115149. 92+452729. 8 


53084 


1368 


503 


8829 


-4.75 


0.04 




SDSS J122545. 87+470613.0 


53117 


1451 


35 


6109 


-5.82 


0.01 


DQ 


SDSS ,1123347.60+125346.1 


53169 


1616 


429 


7152 


-6.36 


0.02 


DQ 


SDSS ,1130945.62+444541.0 


53084 


1375 


391 


8085 


-4.01 


0.03 


DQ 


SDSS J131534. 72+471108. 9 


53062 


1461 


428 


7524 


-5.99 


0.01 


DQ 


SDSS J131930. 66+140137.1 


53112 


1773 


105 


7626 


-5.65 


0.05 


DQ 


SDSS J133127. 04+670419. 5 


51988 


496 


583 


8899 


-4.91 


0.06 


DQ (weak bands) 


SDSS J143144. 83+375011. 9 


53089 


1381 


599 


6173 


-6.97 


0.02 


DQ 


SDSS J152812. 05+513445. 2 


52378 


795 


319 


7531 


-5.69 


0.03 


DQ 


SDSS J153447. 54+414559. 4 


53149 


1679 


616 


7804 


-5.92 


0.02 


DQZ (Call) 


SDSS J161653. 36+392444. 4 


52759 


1336 


572 


7319 


-5.92 


0.02 


DQ 


SDSS J165436. 86+315754. 4 


52791 


1176 


238 


7258 


-5.89 


0.01 


DQ 


SDSS J171341. 76+324009.1 


52413 


976 


623 


7901 


-5.36 


0.01 


DQ 


SDSS J211130. 04-003628. 8 


52431 


985 


35 


7168 


-6.27 


0.05 


DQ 


SDSS J213503. 32+000318. 4 


52468 


989 


198 


6413 


-6.78 


0.03 


DQ 



Here c b and cth are the four observed and theoretical 
SDSS colors, and a the measurement error of the obser- 
vation. We eliminated the majority of objects by using a 
limit of d = 6.5. This limit had been determined from a 
test using the 40 DQs of Fig. ^ The final sample contained 
6952 DQ candidates. 

2.3. Color space for different surface gravities 

We have also fitted the theoretical grids for logg = 7.5 
and 8.5 in the same way. The resulting planes were prac- 
tically indistinguishable from the log<? = 8.0 plane, how- 
ever, the lines of constant temperature or C abundance are 
slightly shifted. This corresponds to the well known fact 
for DQs that the parameters logg and [C/H] are degen- 
erate. They cannot both be determined from photometry 
(nor spectra), because a change of log g can always be 
very nearly compensated with a corresponding change of 
[C/He], and one of the two has to be assumed to be able 
to proceed. Using the fact that most white dwarfs cluster 
around 0.6 M Q , we have assumed logg = 8, as did DB05. 



2.4. Visual classification of the 6952 candidates 

Spectra for all remaining candidates were extracted from 
the SDSS database and inspected visually. More than 95% 
turned out to be quasars, easily recognized by their strong 
and broad emission lines. Of the remaining objects some 
were DA white dwarfs (broad Balmer absorption), some 
galaxies with narrow absorption lines. The DQs, the topic 
of this study, were identified through their molecular Swan 
bands. In total we found 65 DQs, of which 25 were already 
known in the literature and 37 are new detections, and 3 
reclassifications from a spectral type of DC. 

Tables □ and list all 65 DQ with their SDSS names, 
internal identifiers, and atmospheric parameters, deter- 
mined as described in the next section. Most of the objects 
show only carbon features. SDSS J090514. 78+090426.5 
shows a feature near the position of Ha, but nothing 
near the other Balmer lines, and the presence of hydro- 
gen is highly uncertain. SDSS J153447. 54+414559.4 shows 
Call H and K lines with equivalent widths of 2.1 and 
1.0 A. The wings are broad and the distance of the 
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Table 2. 25 rediscovered DQ white dwarfs in the s ame format as Table U 20 of these are analyzed also in DB05. 
Information on the r e maini ng 5 can be found in th e iMcCook fc Sio 3 G22a) catalog rWD 1426+6131. iLiebert et alJ 
(|2003[KlCarollo et alJ {20061. and lHarris etaD l|2003[) . 
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SDSS J000011. 57-085008. 4 


52143 


650 


450 


8042 


-5.46 


0.08 




SDSS J000807. 54-103405. 6 


52141 


651 


199 


7768 


-5.66 


0.06 




SDSS J002531. 50-110800. 9 


52145 


653 


86 


8367 


-4.96 


0.01 




SDSS J015433. 57-004047. 2 


51871 


403 


268 


7435 


-5.89 


0.03 




SDSS J015441. 75+140308.0 


51877 


430 


558 


6511 


-6.89 


0.01 




SDSS J032054. 11-071625. 4 


51924 


460 


236 


6266 


-5.45 






SDSS J033218. 22-003722.1 


51810 


415 


240 


8600 


-4.62 


0.03 




SDSS J090157.92+575135.9 


51924 


483 


600 








Liebert et al. ( 20031 


SDSS J091922. 18+023605.0 


51929 


473 


458 


11566 








SDSS J093537. 00+002422.0 


52314 


476 


461 


4958 


-6.19 


0.02 


Harris et al. C2003) 


SDSS J094004. 64+021022. 6 


52026 


477 


493 


7283 


-5.95 


0.01 




SDSS J095137.60+624348.7 


51943 


487 


227 


8388 


-5.11 


0.10 




SDSS J113359. 94+633113. 2 


52059 


597 


139 


12082 








SDSS J114851. 68— 012612.8 


52056 


329 


578 


9174 


-3.73 


0.01 




SDSS J125359. 61+013925. 6 


52026 


523 


252 


8282 


-4.98 


0.02 




SDSS J123752. 12+415625.8 


53090 


1454 


146 


5846 


-5.51 


0.01 


Carollo et al. (2006') 


SDSS J132858. 20+590851.0 


52411 


959 


504 








Liebert et al. So3) 


SDSS J142728. 30+611026. 4 


52368 


607 


379 


6427 


-6.83 


0.01 


WD 1426+613 


SDSS J144407.25+043446.8 


52026 


587 


418 


9449 


-3.65 


0.05 




SDSS J144808. 07-004755. 9 


51662 


308 


145 


7063 


-6.50 


0.04 




SDSS J154810. 66+562647. 7 


52072 


617 


551 


8119 


-5.46 


0.07 




SDSS J155413. 53+033634. 5 


52023 


595 


373 


6512 


-6.94 


0.03 




SDSS J164328. 54+400204. 3 


52050 


630 


386 


7144 


-6.20 


0.07 




SDSS J165538. 51+372247.1 


52071 


632 


92 


8997 


-4.75 


0.06 




SDSS J205316. 34-070204. 3 


52176 


636 


267 


6382 


-5.45 


0.02 
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Fig. 2. Comparison of effective temperatures determined 
from the photometry with results from DB05 



star is very likely less than 100 pc, implying a photo- 
spheric origin of the lines, probably due to accretion. 
SDSS J094014. 65+090641. 8 sh ows rounded bands similar 
to SDSS J223224.0-074434.3 inlHarris et all J2003I). 



3. Atmospheric analysis of the DQ white dwarfs 

3.1. Photometry 

The parameters T e g and [C /He] were determined by min- 
imizing in the sense of a \ 2 the distance between the ob- 
served point in 4-dimensional color space with model col- 
ors interpolated on the T e g - [C/Hc] grid, \ 2 minimization 
was ob tained with the AMOEBA routine from lPress et alJ 
|I1992Il Since our sample contains 20 objects in common 



with DB05, we compare the derived T e g in Fig. The 
agreement in T e g with those from DB05, which are based 
on photometry and spectra is excellent, confirming that 
photometry alone already gives a reliable T c g. From the 
differences between the two completely independent de- 
terminations we estimate a typical error for T e g of 190 K. 
This is a much more realistic error than the internal er- 
rors from the y 2 routine; it is also very similar to the 
T e fT error of 170 K cited by DB05, and we take this as our 
error estimate for T e g. On the other hand, the derived car- 
bon abundances differ markedly, with a much higher scat- 
ter and also systematic differences. This is not really sur- 
prising, since the colors for different abundances converge 
at the high temperature end and a [C/He] determination 
based on colors alone will be inaccurate to impossible. In 
addition, small uncertainties of the magnitude zero-points 
could shift the whole grid by a few hundreds of a magni- 
tude, resulting in systematic changes. These abundances 
are not used further and therefore not shown here. 



3.2. Spectroscopy 

Fortunately, SDSS provides photometry and spectroscopy. 
The spectra were fitted with the theoretical spectra from 
our standard grid using the same methods used by us 
for other types of white dwarfs I Koester et all l200li 
iHomeier et all I1998|l In principle the spectroscopic fit 
provides T e g and [C/He] (logg = 8 is always assumed). 



However, since we have shown the high reliability of the 
photometric determination, we have kept the temperature 
fixed and determined only the carbon abundance from the 
spectra, leading to smaller errors of these abundances. 
These abundances are given in Table fj] and For the 



D. Koester & S.Knist: New DQ white dwarfs in the SDSS 



5 



i 

o 



-7 -6 -5 -4 -3 

[C/He] 

Fig. 3. comparison of carbon abundances from spectral 
fitting with results from DB05 

five objects with missing entries in the [C/He] column we 
could not obtain reasonable fits within our grid. The spec- 
tral fits for the region around the strongest Swan bands 
of the remaining 60 objects are shown in Figure IA.1I to 
Fig. IA.3I in the appendix; the values are also in Table ^ 
and 

Figure |21 gives a comparison of our final results for the 
carbon abundance with DB05 for the common objects. 
The general agreement is very good, with a small system- 
atic shift of about 0.2 dex, increasing from lower to high 
abundances. The [C/He] abundances of DB05 are on av- 
erage slightly higher than our values. 

4. Results and discussion 

We have identified 40 new DQ white dwarfs and redis- 
covered another 25 from the Data Release 4 of the Sloan 
Digital Sky Survey and presented effective temperatures 
and carbon abundances for 60 of them. This increases very 
significantly the number of DQs with a careful analysis us- 
ing up-to-date model atmospheres. Particularly important 
is the fact that our sample is based on a completely ho- 
mogenous set of photometry and spectra. This is similar 
to the work of DB05, who however included also observa- 
tions from other sources together with 40 SDSS objects in 
their sample. For the objects in common with their anal- 
ysis we find good agreement for the parameters, with a 
slight offset in the abundances of the order of 0.2 dex. 

The major result of this paper is shown in Fig.^J which 
displays the carbon abundance versus effective tempera- 
ture. In addition to our own results (shown as full circles) 
we have here included data from DB05 for those objects 
not included in our sample (open circles) . We clearly con- 
firm the three main conclusions of DB05: 

— The paucity of objects with T e s below 6000 K. Our 
sample adds only one object clearly below the apparent 
cutoff, also carbon-rich as the single cool objects in 
DB05. 
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Fig. 4. Carbon abundance versus effective temperatures. 
Full circles are the results from this paper. We have added 
as open circles the results from DB05 for those objects not 
in our sample 



— The clear sequence for the majority of the objects, from 
low T eff , low [C/He] to high T cff , high [C/He], indicated 
by the lower of the two regression line fits (these fits 
use only the results from this paper). The existence 
of this sequence was first pointed out by DB05 and 
shows the large advantage of using a homogenous sam- 
ple of go od quality observations comp ared to e.g. the 
study of IWeidemann & Koesterl l)l995h . The extrapo- 
lation of the sequence to 10000 K predicts [C/He] = 
-4.0, in agreement with the theoretical calculations of 
iFontaine fc Brassardl {2005) for a He-layer thickness of 
10" 4 of the stellar mass. 

— The existence of a separate population of carbon-rich 
DQs, with an abundance about 1 dex higher than the 
primary sequence. This population seems to form a 
second sequence, which we have tentatively described 
with the upper regression line. 

Based on one single object with a parallax and thus 
mass determination of 1.05 Mq in the upper sequence 
(G47-18 = WD 0856+33), DB05 tentatively suggested 
that all objects of this group are massive white dwarfs 
and might be the descendents of the hot DQ stars, 
which also seem to be massive, alt hough that as s ertion 
is also only based on t wo objects jLiebert et all [2003; 
iMacdonald et all Il998l IWeidemannl |2005|) . Stellar evo- 
lution calculations predict smaller He envelope masses for 
more massive progenitors (IKawai et al., 1988) and might 
lead to the dredge-up occu rring at higher effective temper- 
atures This would also be a natural 
explanation for th e absence of the high-mass tail in the DB 
mass distribution ijBeauchamp et all Il996t iLiebert et all 
120031) . 

While these conclusions seem rather speculative at 
present, it is obvious that further study of the large num- 
ber of DQs with excellent observations coming from SDSS 
and other large scale survey will provide important clues 
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for the remaining open question of the origin and evo- 
lution of white dwarf surface compositions and spectral 
types. Of particular importance are parallax determina- 
tions for DQs on the high abundance sequence. 
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Fig. A.l. Spectral fit for DQ white dwarfs. Thin (blue): observed; thick (red): model. Vertical axis is relative intensity 
on an arbitrary scale. 
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Fig. A. 2. Spectral fit for DQ white dwarfs; continued from Fig. IA.ll 
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Fig. A. 3. Spectral fit for DQ white dwarfs; continued from Fig. IA.2I 



